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ABSTRACT 



The dynamics of a proposed marine nuclear propulsion plant, 
under the condition of sudden changes in load, are studied 
by both an analog and digital simulation. Two inherent reac- 
tor characteristics, internal reactivity feedback and the 
existance of delayed neutron emitters, which effect power 
plant stability and controllability are considered in detail. 

Comparison of the simulation results with practical data 
indicates that the simulation represented the dynamics of a 
marine nuclear power plant with a sufficient degree of accuracy 
to be extremely useful in the study of marine nuclear power 
plants . 

Evaluation of the simulation results indicate that the 
proposed power plant is inherently stable when operating in 
the normal region. The desirability of control schemes are 
discussed and a control scheme utilizing a constant average 
coolant temperature program is implemented. This external 
control scheme significantly improved the response of the 
system. 



2 



TABLE OF CONTENTS 



I . INTRODUCTION 

II. MATHEMATICAL MODEL 

A. REACTOR KINETICS 

B. HEAT TRANSFER FROM FUEL AT T„ TO COOLANT AT T_ T7 - 

r AV 

C. HEAT BALANCE FOR COOLANT IN THE CORE 

D. HEAT TRANSPORT DELAYS 

1. Time Delay from Reactor to Boiler 

2. Time Delay from Boiler to Reactor 

E. HEAT BALANCE FOR HEAT EXCHANGER (BOILER) 

1. Primary-Side 

2. Secondary-Side 

F. FINAL MODEL WITH SPECIFIC SYSTEM PARAMETERS 

1. Reactor Kinetics 

2. Heat Transfer from Fuel to Coolant 

3. Heat Balance for Coolant 

4. Time Delays 

5. Heat Balance for Steam Generator Secondary- 
Side — 

III. INHERENT REACTOR CHARACTERISTICS 

A. PROMPT AND DELAYED NEUTRONS 

1. Kinetic Equations Assuming No Delayed 

Neutrons 

2. Reactor Kinetic Equations with a One Group 

Delayed Neutron Emitter 

3. Comparison of Results 



12 

15 

15 

18 

21 

22 

22 

23 

23 

24 

25 

26 

26 

26 

26 

27 

27 

28 

28 

29 

30 

30 



3 



B. INTERNAL REACTIVITY FEEDBACK 

1. Temperature Effects 

2. Fission-Product Poisoning 

3. Net Internal Reactivity 

C. TOTAL EFFECT OF DELAYED NEUTRON AND NEGATIVE 

REACTIVITY FEEDBACK 

IV. DIGITAL AND ANALOG SIMULATION OF POWER PLANT WITHOUT 

EXTERNAL FEEDBACK 

A. DIGITAL SIMULATION 

1. Load Changes 

2. Results 

3. Discussion of Results 

B. ANALOG SOLUTION 

1. Prompt- Jump Approximation 

2. System Variables to Computer Variables 

a. Reactor Kinetics 

b. Heat Transfer from Fuel to Coolant 

c. Heat Balance for Coolant 

d. Time Delays 

e. Heat Balance for Steam Generator 

Secondary-Side 

3. Comparison of Analog and Digital Solution 

V. CONTROL 

A. STEADY STATE PROGRAMMING 

1. Constant T. tt 

AV 

2 . Constant T g 

3. Comparison of Constant T and Constant T 

Programs 

B. DESIGN OF EXTERNAL CONTROL SYSTEM 

1. Reactivity Rate 



33 

34 

34 

35 

35 

38 

38 

38 

38 

39 

39 

39 

40 

40 

41 

41 

41 

42 

42 

46 

46 

47 

47 

48 

49 

49 



4 



2. Automatic Controller 



50 



3. Evaluation of External Controller 51 

VI. CONCLUSIONS AND RECOMMENDATIONS 56 

A. CONCLUSIONS 56 

B. RECOMMENDATIONS 57 

APPENDIX A SPECIFIC SYSTEM EQUATIONS 59 

APPENDIX B TIME RESPONSE OF POWER PLANT 62 

APPENDIX C PREPARATION OF SYSTEM EQUATIONS FOR USE 

ON ANALOG COMPUTER 6 7 

APPENDIX D CONTROLLER SIMULATION 70 

BIBLIOGRAPHY 77 

INITIAL DISTRIBUTION LIST 79 

FORM DD 147 3 80 



5 



LIST OF TABLES 



4-1 COMPARISON OF ANALOG AND DIGITAL SOLUTION FOR A 

STEP LOAD CHANGE OF 20 PERCENT 4 4 

4- 2 COMPARISON OF ANALOG AND DIGITAL SOLUTION FOR A 

STEP LOAD CHANGE OF 40 PERCENT 45 

5- 1 COMPARISON OF SIMULATED POWER PLANT WITH AND WITHOUT 

EXTERNAL CONTROLLER FOR A STEP LOAD CHANGE OF 20 
PERCENT 53 

5-2 COMPARISON OF SIMULATED POWER PLANT WITH AND WITHOUT 
EXTERNAL CONTROLLER FOR A STEP LOAD CHANGE OF 40 
PERCENT 54 



6 



LIST OF DRAWINGS 



1- 1 BLOCK DIAGRAM OF BASIC NUCLEAR POWER PLANT 13 

2- 1 TEMPERATURE PROFILE OF SINGLE FUEL ELEMENT 18 

2-2 TEMPERATURE PROFILE ACROSS FUEL ELEMENT AND 

CLADDING 20 

2- 3 TEMPERATURE PROFILE FOR COUNTER FLOW HEAT EXCHANGE - 24 

3- 1 TIME RESPONSE OF REACTOR CONSIDERING ONLY PROMPT 

NEUTRONS 31 

3-2 TIME RESPONSE OF REACTOR WITH PROMPT AND DELAYED 

NEUTRONS 32 

3- 3 TIME RESPONSE OF REACTOR WITH DELAYED NEUTRONS AND 

REACTIVITY FEEDBACK 36 

4- 1 ANALOG COMPUTED DIAGRAM FOR BASIC POWER PLANT 43 

5- 1 CONSTANT T^ v PROGRAM 47 

5-2 CONSTANT T g PROGRAM 47 

5-3 BLOCK DIAGRAM OF POWER PLANT WITH EXTERNAL 

CONTROLLER 52 

B-l TIME RESPONSE OF SIMULATED POWER PLANT TO A 40 

PERCENT STEP LOAD CHANGE (DIGITAL) 63 

B-2 TIME RESPONSE OF SIMULATED POWER PLANT TO A 40 

PERCENT STEP LOAD CHANGE (ANALOG) 65 

D— 1 BLOCK DIAGRAM OF CONTROL ROD POSITIONING SYSTEM 71 

D-2 TIME RESPONSE OF CONTROL ROD POSITIONING SYSTEM 73 

D-3 TIME RESPONSE OF CONTROL ROD POSITIONING SYSTEM 74 

D-4 TIME RESPONSE OF POWER PLANT WITH EXTERNAL 

CONTROLLER 75 



7 



TABLE OF SYMBOLS AND ABBREVIATIONS 



SYMBOL 

A 

A 1 



c 

c (0) 



c c 

c f 



E 

k 

k l 

k 2 

K 

K R 

k m 

k t 



l 

ra c 

III, 



m 



DEFINITION 

Area 

Constant of proportionality 
relating power and neutron 
density. 

Density of delay neutron emitters 

Density of delayed neutron 
emitters at time t = 0 

4 

Specific heat of the coolant 

Specific heat of the fuel 

Specific heat of the steam 

Normalized density of delayed 
neutron emitters . 

Arbitrary constant 

Excess reactivity 

Coefficient of heat transfer 

Coefficient of heat transfer 

Excess reactivity in dollar 
units 

Reactivity associated with rod 
position 

Motor constant 

Throttle constant 

Prompt neutron lifetime 

Mass flow rate of the coolant 

Mass of the fuel 

Mass flow rate of the steam 



8 



SYMBOL 


DEFINITION 


n 


Neutron density 


n (0) 


Neutron density at time t = 0 


N 


Normalized neutron density 


P 


Power 


q 


Heat transferred 




Heat transferred to the coolant 
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I . INTRODUCTION 






Nuclear reactor power plants for use in marine application 
differ from stationary power plants in that they must respond 
quickly to power demands which may vary throughout a wide 
range in a relatively short period of time. There are maneu- 
vering situations where the shipboard power plant may be 
required to respond to load changes of up to 50 percent in a 
few seconds. In addition to the operational differences, the 
mobile nature of the nuclear power ship coupled with the size 
and weight limitations which must be imposed upon the reactor 
shielding require that safety be a primary consideration in 
design and control of the marine power plant. 

There are several references [Harrer 1963 and Glasstone 
and M. C. Edlund 1952] that treat the nuclear reactor and 
nuclear reactor control from the standpoint of the reactor 
operating at constant output power. Schultz [Ref. 3] discusses 
the problems of reactor control from the standpoint of a 
reactor operating with varying power demands as well as when 
operating under steady power out conditions. King -4-Reir— - 4~Y 
considers the problems associated with nuclear power systems 
largely from the position of a mechanical engineer. 

This thesis investigates a marine nuclear power plant and 
is intended as an initial look into the dynamics of this type 
of propulsion system. The system dynamic equations are 
'developed and the resulting fundamental balance equations are 
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simplified by invoking reasonable first approximations. 
Initially a power plant without an external controller is 
studied. After the initial investigation a control scheme is 
formulated and implemented. The plant response with controller 
is compared with the plant without controller. 

The power plant, represented in block diagram form in 
Figure (1-1) , utilizes a pressurized water reactor fueled with 
low enriched uranium dioxide, heat exchanger for the conver- 
sion of reactor power to steam and a turbine which is repre- 
sented by the changing load. The reactor and one side of the 



I 




I 

I 



Fig. (1-1) . BLOCK DIAGRAM OF BASIC NUCLEAR 
POWER PLANT 



heat exchanger represent the primary loop while the other side 
of the heat exchanger and turbine represent the secondary loop 
in the system. 
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The system, although hypothetical, is similar to the 
plant on board N. S. Savannah. Many of the initial (before 
load changes) temperatures and coefficients closely approxi- 
mate those of Savannah — 6 , 7 ~]~. In order to show the valid- 
ity of the mathematical model used in this study, a qualitative 
comparison between the response of the simulation and both 
those of Savannah and a proposed marine power plant is used 
[ 8 ). 

The organization of this thesis follows the format as 
described below: Chapter (2) contains a development of- the 

mathematical model used for the evaluation of system perfor- 
mance; Chapter (3) contains a general discussion of two 
inherent reactor characteristics, delayed neutrons and internal 
reactivity feedback, that affect the stability of the power 
plant; Chapter (4) presents an analysis of the analog and 
digital computer simulations of the power plant without 
external control and also contains a comparison of the results 
of the two simulations; Chapter (5) discusses some general 
aspects of nuclear power plant control and shows the specific 
control scheme used in this evaluation; Chapter (6) is a 
discussion of the conclusions arrived at by this study and 
also discusses areas of further interest. Appendices A 
through D supplement these chapters. 
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II. MATHEMATICAL MODEL 



The mathematical model for the plant under consideration 
is developed by obtaining the model for each system component 
individually. The individual models together with the values 
of the specific system parameters are summarized at the end of 
this chapter. During the course of the development, it is 
necessary to make several simplifying assumptions. These 
assumptions introduce some inaccuracies in the system model. 
"However, the errors involved in these assumptions are 
usually less than the amount of uncertainty in the engineering 
value of the coefficients used." [3] The generalized equa- 
tions given below are listed in [Ref. 9] . 



A. REACTOR KINETICS 



The point model kinetic equations of a chain reacting 
nuclear pile have been developed in the literature [10, 11] 
and shall be used for the starting point in the development 
of the plant model to be studied. 




6k-B 
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6 
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A i C i 



( 2 - 1 ) 



c. 
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h 



( 2 - 2 ) 



Equations (2-1) and (2-2) are state equations describing a 
nuclear reactor operating at some steady state power condition 
The variables n, c, and 5k represent neutron density, density 
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of delayed neutron emitters and excess reactivity. The other 
terms are constant coefficients described ea~-pa.gfe~?. 

The indexed terms in equation (2-1) and (2-2) indicate a 
multi-group representation of the delayed neutron emitters. 
This six group representation is valid for a detailed study 
of reactor dynamics, however, in the scope of this study the 
reactor is considered as a power source which is investigated 
from the terminal point of view. Therefore, without loss of 
generality, and for the convenience and simplicity (i.e., 
equipment limitation in the analog system) , the kinetic 
equations are represented by a one group delayed neutron 
approximation. Equations (2-1) and (2-2) thus become: 



The magnitudes of the variables n and c are far greater 
than the magnitudes of the other system variables. in -order 
to facilitate the use of the analog computer, equation (2-3) 
and (2-4) will be normalized by dividing through by their 
steady state values n(0) and c(0). 




(2-3) 




(2-4) 
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n (0) 
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(2-5) 
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c (0 ) 



Ac 



( 2 - 6 ) 



16 



Next defining: 



N 



n 

n C o) 




and K 



6k 

B 



noting that in steady state: 



N = D = 1 , N = D = K= 0; and n (0 ) = X c (0) 
Equations (2-5) and (2-6) thus become: 

N = |- (KN - N + D) (2-7) 

D = A(N - D) (2-8) 



The term K above is the effective net reactivity existing 
within the reactor. In steady state there is a constant 
neutron population and therefore K must be zero. To establish 
K at this zero value the reactivity attributed to rod position 
must be off set by the inherent internal negative feedback 
(internal feedback is discussed in some detail in Chapter 3) . 
Thus the algebraic equation for K is: 



K = k r + 



“f t f + 



a T mt 

c AV 



(2-9) 



Tp is the average fuel temperature and T ^ is the average 
coolant temperature. The coefficients and are coeffi- 
cients of reactivity feedback. It should be noted that a p 
and a c are negative .numbers . 
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B. HEAT TRANSFER EQUATIONS FROM FUEL AT T TO COOLANT 
AT T av 

For a given change in temperature the fuel will absorb 
heat in the amount of 

q f = m f c f CT 2 " T l } T 2 > T 1 (2-10) 

Where is the heat absorbed by the fuel, m^ is the mass of 
the fuel, is the specific heat of the fuel and T£ - T^ is 
the change in fuel temperature. The time rate of change of 
the heat absorbed in this manner is: 

q f = m f c f (t 2 " T i } (2_11) 

Consider one reactor fuel element of Figure (2-1) . The 
temperature profile of the fuel varies along the axial length 
Z and depends upon : 

1) the heat input from the reactor. 

2) the coolant flow rate, which is constant. 




Fig. (2-1). 



TEMPERATURE PROFILE OF 
SINGLE FUEL ELEMENT. 



The first of several assumptions must now be made. Assume 
that the temperature profile does not change shape or form 
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during transient conditions and that in the small signal 
analysis and m^ are constant. Then at any axial position 
Z the rate of heat absorption by the fuel is dependent upon 
the time derivative of fuel temperature. Further considering 
changes from steady state where T^ is a constant, it then 
follows : 



= 



( 2 - 12 ) 



The general heat transfer equation is [12] 



q = k-i V' 



(2-13) 



axial (z direction) flow can be neglected in comparison to 
radial flow. For a homogeneous right cylinder fuel element 
equation (2-13) becomes : 
k. 



• K 1 2 / 2T \ 

q = — 27 ( r 2F ) 



(2-14) 



employing the power balance for each fuel element where Power 
= time rate of change of heat 



Power generated 
by the fission 
process 



Power 

absorbed 1 + 
t by the fuel, 



Power 

transferred 
to the coolant 



Considering the case where k^ is independent of space 
variables, the (differential) power balance equation across a 
differential radius r is: 



^ m f C f ^F ^1 




1 2T 
r 2r 



(2-15) 



19 



or rearranging 



T = l - [ 2 2 T 1 2T\ 

f m f c f m f c f \ 2t 2 r 2r ) 



(2-16) 



Since and k^ are both, functions of temperature, the above 
equation is a non-linear partial differential equation. 
Following the assumption that the temperature profile does 
not shape or form and that 2T/2r is a constant as shown in 
Figure (2-2) 
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Ar 



t f t av 
r 2" r l 



2 

2 T 
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2r 



= 0 



TAV = 



T„+T 
H c 



Fig. (2-2). TEMPERATURE PROFILE ACROSS 
FUEL ELEMENT AND CLADDING. 



thus equation (2-16) becomes 

k. 



T = — 
F m 



f C f 



m f c f 



t f t av \ 
r 2 " r l 



(2-17) 



The power generated is a function of the neutron population 
[1] . Thus : 



P 



nr 



f (N) 



A 1 N 
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The term 



has the units of sec 



-1 



and can 



m f c f r(r 2 ' r l> 



be related to a time constant t^. Therefore equation (2-17) 
becomes : 



T f = A 1 N 




T ) 
AV ; 



(2-18) 



C. HEAT BALANCE FOR COOLANT IN THE CORE 

Under the same assumptions as stated above the rate of heat 
absorbed by the coolant as it passes through the fuel elements 
is : 

- m c h (T H - V < 2 - 19 > 



Where a is the rate of heat transferred to the coolant, m 
'C c 

is the mass flow rate of the coolant, c is the specific heat 
of the coolant, T is the hot leg temperature and T., is the 
cold leg temperature. Employing a power balance as above 

( Power transferred ' _ Power absorbed \ + / Power transferred 
to the coolant / x by the coolant J \ to the boiler 




r 






T ) 
C 1 



( 2 - 20 ) 



Where t is the time for a unit volume of coolant to flow 
o 



through the reactor. 



'AV 



k i ct av - V 

rCr 2 ' r l> m c C c 



- T- (t h - T c> 

o 



( 2 - 21 ) 
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can be related to a time constant x 



again 
Thus : 



r(r_ - r, ) rn c 
Z ICC 



*av - tT < t av - t f» - r- (t h - T c> 

2 o 



Rearranging 



T T^. 2 T 

T - F AV c 

AV T 0 T _ 1 T 

Z Z O 



where 



x x _ 
o 2 

x + 2x 0 
o 2 



2 * 



( 2 - 22 ) 



(2-23) 



D. HEAT TRANSPORT DELAYS 

There are time delays associated with the transfer of heat 
from one portion of the plant to another. In developing the 
model for the time delays the following assumptions must be 
made : 

1) there is no heat loss in the piping between the reactor 
and the boiler. 

2) there is no mixing of fluids in the system piping. 

1 . Time Delay From Reactor to Boiler 

The boiler inlet temperature , T , is a delayed version 

DX 

of the reactor outlet temperature Ty. 

t bi - t h (t - t 3> (2 ‘ 24) 
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Where T is the boiler inlet temperature and x _ is the 
riJL J 

transport time between reactor and boiler. Rearraning 

T H = T BI (t + T 3> (2 ' 25) 

Expanding in a Taylor Series 

T H = T BI + 1 3 ^BI + * • * (2-26) 

Solving for T BI 

*BI " 7J CT BI - 2 T AV + T C> < 2 - 27 > 

2 . Time Delay From Boiler to Reactor 
Following the above reasoning 

T c - 77 <T c - 2 t b + t bi» (2 - 28) 

Where T g is the average temperature in the boiler. 

E. HEAT BALANCE FOR HEAT EXCHANGER (BOILER) 

The development of the model of a counter flow heat 
exchanger in which there is a change of phase (water to steam) 
in one flow path is a difficult task. It is not the intention 
of this chapter to present a rigourous derivation of the equa- 
tions used in the system model, but rather it is to give 
insight into the physical relationships from which these 
equations were obtained. Therefore a more heuristic argument 
will be used in the development of the steam generator model. 
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TEMPERATURE 




Fig. (2-3) . TEMPERATURE PROFILE FOR COUNTER 
FLOW HEAT EXCHANGER. 

There are several additional assumptions which must now be 
made : 

1) There is no heat lost or gained with the surroundings. 

2) The specific heat of the individual fluids remains 
constant . 

3) The overall coefficient of heat transfer (k^) 
uniform through the heat exchanges. 

4) From Figure (2-3) the mean temperature difference 
throughout the heat exchanges T m = t b - T S 

1 . Primary-Side 

Consider a small section of heat exchanger having a 
transfer area dA. Then the incremental heat transfer is: 

dQ » k_ dA AT = c m dT (2-29) 

2 n c c c 

integrating yields 

Q = k 2 A[T b - T s ] = c c m c [T bi - T b0 ] (2-30) 



24 



Since k~, A, c , and m are all constant, equation (2-30) 

" C 

reduces to 



Where T g is the steam temperature and A is the area of the 
heat exchanger. 

2 . Secondary-Side 

Energy balance in th.e secondary-side consist of the 
amount of power transferred to the secondary, k^A(T B - Tg) , 
equal to the power absorbed within the boiler plus the power 
which is transferred to the turbine. To bring this relation 
into a linear, ordinary differential equation some additional 
assumptions are required. These are: 



1) The thermal capacity of the boiler metal, steam 
and coolant can be lumped. 

2) The power demand on the boiler can be measured as 
a fraction of normal load. ip will represent this 
condition with tJ; = 1 representing the normal load. 

3) The heat exchanged in the boiler occurs at a single 
point within the heat exchanger. 



Therefore the flow of heat delivered to the turbine may be 
represented as the product of the throttle opening, represented 
hy K t , and the steam pressure. Since the steam pressure is 



to represent steam pressure. Thus the heat balance for the 
secondary side becomes: 



2(T b - T s ) = E(T bi - T bq ) 



(2-31) 



2P, 



upon the steam temperature, the term T can be used 
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Where rn and c are the mass flow rate of the steam and the 
s s 

specific heat of the steam respectfully. Rearranging 



k.. A 

• I 

T = = FT 

S m c + m c 1 B 



T s> + TS 



k t * 



c + m c ) 
c c s s 



2P 

S 

— T 

2T„ S 



s 

(2-33) 



Which reduces to 



1 

= - - C-T b + T s + K* ^ 2rf 



(2-34) 



F. FINAL MODEL WITH SPECIFIC SYSTEM PARAMETERS 

The development of the specific equations for the system 
is carried out in detail in Appendix A, they are summarized 
here for convenience. 



Reactor 


Kinetics 




N = 6.4 


(KN-N+D) 


(2-35) 


D = .1 


(N-D) 


(3-36) 


* 

1! 

* 


- 1.57 x 10 -2 T p - 3.14 x 10“ 2 T av 


(3-37) 



where N(0)=1, D(0)=1 and K (0) =24 . 7275 . 

K 

2 . Heat Transfer From Fuel to Coolant 

T p = - | (T p - T av ) + 30.0 N (2-38) 

where T^ (0) =565°F. 

r 

3 . Heat Balance for Coolant 

T av = 10 T f - 50 T av + 40 T c (2-39) 

where T AV (0) =505°F. 
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4 . Time Delays 



T 




(2-40) 



where T BI (0)=520°F 



T 



c 




(2-41) 



where T c (0)=490°F. 

5 . Heat Balance for Steam Generator Secondary-Side 



T 



S 



j (- | T bi + (0.667 + 0.0665 ij>) T g ) 



(2-42) 



where Tg(0)=565°F. 

G. METHOD OF SOLUTION 

Because of the multiplication of the time dependent 
variables K and N in equation (2-7) and because K is a function 
of N the entire set of dependent equations which describe the 
plant dynamics is non-linear. Since (in general) non-linear 
differential equations do not lend themselves readily to 
analytical solutions, computer simulation will be used to 
solve these equations. 
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III. INHERENT REACTOR CHARACTERISTICS 



There are several important characteristics of nuclear 
power plants which greatly effect the stability and control- 
lability of the system. The two most significant factors 
relating to the system performance are the existance of 
delayed neutron emmiters and the negative internal reactivity 
feedback associated with temperature changes within the 
system. Both of these phenomena are discussed below in some 
detail. 

A. PROMPT AND DELAYED NEUTRONS 

At the instant of fission there are neutrons (prompt) 
which are produced instantaneously. Fortunately, for control 
purposes, a small portion (g) of all neutrons produced during 
the fission process are delayed in time [1, 2]. These delayed 
neutrons are the result of the decay of fission products. As 
noted previously these delayed neutrons are grouped in six 
distinct groups with different delay constants (X) . Analysis 
of the effects which delayed neutrons exhibit on the reactor 
kinetics can best be shown by the evaluation of the kinetics 
equations (2-3) and (2-4) with and without the delayed 
neutrons considered. 

The analysis of the prompt and delayed neutron response 
which follows is a study of the reactor operating at a power 
level low enough so that feedback effects (discussed later 
in this chapter) are negligible. An analagous situation in 
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a power plant would be In the start-up power range of the 
reactor . 

In this evaluation the input to the reactor is in the 
form of a perturbation to the steady state reactivity k. The 
input was limited to a maximum value (6) of (0.0064) , for 
when 6k exceeds 3 the reactor is in a response condition 
described as "prompt critical." "Prompt critical" refers to 
the fact that the nuclear fission chain reaction can be main- 
tained by means of prompt neutrons alone. If this condition 
occurs neutron density, and hence power increase rapidly 
from the instant of input, making the reactor difficult to 
control and may cause what is referred to as a start up 
accident. The "prompt critcial" operation is avoided in 
practice. 

1 . Kinetic Equations Assuming No Delayed Neutrons 

Under the assumption of no delayed neutrons, equations 
(2-3) and (2-4) , which describe the reactor kinetics in steady 
state with a one group delayed neutron approximation, reduce 
to the single equation: 




Solving for n 



n = n (0) e 



6k 

l 



(3-2) 



Dividing by n(0), substituting the system parameters found in 
Appendix A, and assuming 6k = 0.003 equation (3-2) now becomes 

N = e 3t (3-3) 
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2 . Reactor Kinetic Equations with a One Group Delayed 
Neutron Emitter 

Under the one group delayed emitter assumption, 
equations (2-3) and (2-4) are unchanged. With the substitution 
of system parameters from Appendix A and from the solution of 
these equations carried out in Ref. [11], the solution of 
these equations for N is : 



N 



1.88 e +0 ’°884 t 



0.88 



e 



-3.4 



t 



(3-4) 



3 . Comparison of Results 

The effect of the delayed neutrons is readily seen 
by the comparison of equations (3-3) and (3-4) . With no 
delayed neutrons the reactor neutron population has increased 
by a factor of 8,100 in 3 seconds, however when delayed 
neutrons are present the reactor neutron population has only 
increased by a factor of 2.5 in 3 seconds. Plots of N versus 
time for both of these situations are shown in Figures (3-1) 
and (3-2). As time passes (i.e., time greater than 10 seconds) 
the first term of equation (3-4) becomes dominant and the 
response will follow a positive exponential increase. However, 
in the time span of interest (a small time after the distur- 
bance) the affect of the delayed neutrons is to greatly reduce 
the rate of change of neutron flux. It is the natural occur- 
ance of these delayed neutrons which makes it possible to 
operate a nuclear reactor with the degree of safety required. 
Even as small a fraction as 0.64 percent of the total neutron 
population acting as delayed neutrons has the effect of making 
the entire reactor control problem greatly simplified [2, 3]. 
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It is considered pertinent to re-emphasize that the 
above .discussion of prompt and delayed neutrons was conducted 
for a reactor in the startup range. The response of a 
reactor operating in the normal power range (where reactivity 
feedback is strong) is significantly different than that of 
the reactor operating at low power. This difference will be 
demonstrated later in this chapter. 

It should be pointed out however, that the delayed 
neutron phenomenon also has an undesirable facet. While it is 
readily obvious that it restricts the growth of the neutron 
population for positive 6k inputs, it also restrains the 
reduction rate of neutron population which would result from 
a negative 6k input as the power demand is reduced during 
normal operation or during shut down. 

B. INTERNAL REACTIVITY FEEDBACK 

In the operation of a nuclear reactor the power generated 
is directly proportional to the neutron population [1]. In 
order for the neutron population to remain constant at some 
steady state power level, equation (2-7) must be equal to 
zero. For this situation to occur, K, the net internal 
reactivity which is the sum of both the reactivity attributed 
to the control rod position and the reactivity due to internal 
feedback, must be zero. Changes in the feedback reactivity 
depend upon temperature variations and fission-product — 
poisoning which are discussed below 

To insure the stability of the reactor the overall coeffi- 
cient of reactivity feedback must be negative [13] . 
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1. Temperature Effects 



When there is a power change in a nuclear reactor 
there is considerable additional energy released. Since it 
is not possible to transfer all of the liberated energy to 
the load , there are temperature fluctuations throughout the 
reactor and associated coolant. 

Temperature variations affect reactivity in two 
distant ways: first, the mean energy of the thermal neutrons 
and hence their nuclear cross-section vary with temperature; 
and, second, the mean free path length and leakage probabilities 
are functions of temperature [2, 3]. These effects result in 
a net negative reactivity. A reactor power plant of the type 
described above, having a negative temperature coefficient of 
reactivity will tend to be self-controlling in response to 
load changes and limited external reactivity variations [14] . 

In hetrogeneous reactors (the type under study here) 
the fuel and the coolant are at different temperatures and at 
different physical states. Therefore, there will be specific 
effects associated with both the fuel and coolant [15] . 

Temperature effects upon reactivity manifest themselves 
in temperature coefficients of reactivity. (a p and a c ) . 

2 . Fission-Product Poisoning 

When a reactor is operating at power certain fission 
products, notably Xenon 135 and Samarium 149, are accumulated. 
These products have large nuclear cross sections and absorb 
the excess neutrons, which are vital in maintaining the chain 
reaction. The rate of formation of these products is 
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dependent upon the rate of fission occuring at the specific 
power level. Eventually the rate of formation and the loss 
of the absorbing nuclei become equal and equilibrium exists. 
However, when the reactor is shut down, these substances, 
which result mainly from the decay of other fission products, 
continue to grow and may reach their highest concentration 
several hours after shut down occurs. Since the production 
of these fission products continues after shut down, it is 
entirely possible that there may be enough negative reactivity 
remaining to restrict the power build up of the reactor until 
these products decay to a level where the inherent positive 
reactivity in the control rods can overcome their adverse 
reactivity effects. 

The above fission ^product phenomenon does not essent- 
ially affect the performance of a reactor operating at power. 
For this reason these effects are not considered in the scope 
of this study. 

3 . Net Internal Reactivity 

The net reactivity in dollar units, due to both rod 
position and to temperature can be expressed as: 



K - k r + 



“f t f + 



a T 
c AV 



C. TOTAL EFFECT OF DELAYED NEUTRONS AND NEGATIVE REACTIVITY 
FEEDBACK 



The total effect of the delayed neutrons and negative 
reactivity can best be illustrated in Figure (3-3) . These 
plots are for the specific plant under study. The power plant 
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3b 



disturbance is a step input of magnitude K = 0.469 dollars 
(i.e. , <5k = 0 .003) . 

The input causes the reactor's neutron density to rise 
sharply from a steady state unit value to a value of 1.6. If 
there were no delayed neutrons this neutron density would 
rise faster and obtain a larger value. The negative reactivity 
feedback acts rapidly and limits the excursion of the net 
reactivity to a maximum value of 0.38. After approximately 
60 seconds, the reactivity feedback has adjusted the internal 
reactivity to zero. The steady state effect of the disturbance 
is to raise the neutron density to a level of 1.02. 
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IV. DIGITAL AND ANALOG SIMULATION OF POWER PLANT 
WITHOUT EXTERNAL CONTROL 



A. DIGITAL SIMULATION 

The system equations were programmed and solved on the IBM 
360-67, using the Naval Postgraduate School Library Program 
(DRKGS) for the integration. This integration method utilizes 
a forth-order "Runge-Kutta" solution with double precision 
arithmetic . 

Time response of the system to load changes of +40 percent 
are shown in Appendix B and are tabulated for all experimental 
load changes at the end of this chapter. 

1. Load Changes 

The system was subjected to standardized step load 
changes disturbances of + 20 percent and + 40 percent. These 
load changes were selected because a 20 percent load change 
represents a change which a maritime power plant should be 
able to respond to quickly, and a 40 percent load change 
represents a feasible change which could be encountered during 
an accidental situation. 

2 . Results 

The results of the simulation of this hypothetical 
system were compared against the time response of marine power 
plant designs of similar size [6, 8]. It was found that the 
simulation used in this study compared well with the referenced 
systems. The steady state (after the load distrubance) values 
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of the system variables coincided well with their counterparts. 
The simulated system demonstrated a slightly more oscillatory 
behavior than the actual and proposed systems referenced. 

This oscillatory behavior is attributed to the simplified 
system model used in this study. Even with the small oscil- 
latory differences, the simulation represents the actual 
system with a sufficient degree of accuracy that it is quite 
suitable for an initial study of marine reactor power plants. 

3 . Discussion of Results 

a) The power plant, without benefit of an external 
control system, was shown to be stable to all subjected load 
changes. The system displayed a second order well damped 
response to load changes . 

b) Temperature changes, including overshoots, are not 
severe and should not cause thermal stress problems. 

c) System response, although sluggish, is fast enough 
for use as a maritime power plant. 

B. ANALOG SOLUTION 

The analog simulation of the system was carried out on 
the Comcor CI-5000 computer. This computer has a + 100 volt 
dynamic range. The CI-5000 computer is interfaced with the 
XDS-9300 digital computer. This hybrid pair will be used 
later in the control section of this study. 

1 . Prompt-jump Approximation 

Figure (B-l) of Appendix B shows that the time rate 
of change of N is not large. For ease of analog simulation 
(i.e., to reduce the number of high gain amplifiers required 
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in the simulation) equation (2-7) can be modified as follows: 



= KN - N + D 

6 . 4 



(4-1) 



making the approximation that N = 0 and solving for N 




(4-2) 



This approximation is called the prompt- jump approximation 
since any change in K is reflected instantaneously in N. The 
system equations using this assumption were simulated on the 
IBM 360-67 and it was found that this approximation does not 
introduce errors greater than 2 percent in the system response. 

2 . System Variables to Computer Variables 

The system variables such as temperature and relative 
neutron population must be related to analog computer variables 
which are voltages . The system equations and initial condi- 
tions in terms of computer variables are developed in Appendix 
C. These equations are summarized below for convenience: 
a. Reactor Kinetics 




(4-3) 



D = 0.1 (N-D) 



(4-4) 



K = 3.5 K 



R 



0.471 T 



F 



0.942 T 



AV 



(4-5) 
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where 



N(0) = 10 volts 

5(0) = 10 volts 

K (0) = 35.36 volts 

K 

b. Heat Transfer From Fuel to Coolant 

T F = - 0.5 [T p - T av ] + 0.5 N (4-6) 

Where 

T„ (0 ) = 94.17 volts 

r 

c. Heat Balance for Coolant 

* AV = -50 T av + 10 T p + 40 T c (4-7) 

Where 

T av ( 0) = 84.17 volts (4-8) 

d. Time Delays 

*BI = -°' 5 ' ? BI - 2 f AV + V (4 ' 9) 

T c = °.S [-f c + |T S + |t bi ) (4-10) 

Where 

T nT (0 ) = 86.67 volts 

-b JL 

9 

T (0) = 81.67 volts 

c 
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e. Heat Balance for Steam Generator Secondary-Side 



T s = 10.133 T - 0.1324 T g - 0.0027 ip T g ] (4-11) 
Where 

T CO) = 79.17 volts 

The analog computer diagram for the solution of these equations 
is shown in Figure (4-1) . 

3 . Comparison of Analog and Digital Solution 

Below is shown two tables comparing the results of the 
digital and analog computer solutions . Because of the inherent 
resolution quality characteristics of the digital computer, 
the digital solution will be used as the base for determining 
errors. Settling time is defined as the amount of time 
required for the magnitude of K(t), the net internal reactivity 
to return to a value of + 0.005 after the applied step load 
disturbances . 

From the tabulated results it is apparent that the 
analog and digital solutions are in such close agreement that 
the analog solution is adaquate for the remainder of this study 
Henceforth the digital computer is used as an auxiliary compu- 
tational tool to solve specific problems but is not used in 
the evaluation of the overall simulation. 

The value of the system variables obtained from the 
analog simulation were obtained through Analog to Digital 
conversion and printed out from the XDS-9300 digital computer. 
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-50 




Fig. (4-1) . ANALOG COMPUTER DIAGRAM FOR BASIC POWER 
PLANT. 
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A ip = + 0 . 2 


t> 

■€- 

II 

1 

o 

• 

to 


ANALOG 


DIGITAL 


% ERROR 


ANALOG 


DIGITAL 


% ERROR 


K (MAX /MIN ) 
$ 


0.068 


0.066 


3 


-0.081 


-0.084 


3.58 


N ( MAX/MI N) 


1.21 


1.20 


0.83 


0.81 


0 . 79 


2.53 


N (SS) 


1.21 


1.19 


1.68 


0.82 


0.80 


2 . 5 


T (MAX/ 

MIN) 
° F 


573.0 


572.7 


0.052 


557.1 


556.8 


0.054 


T (SS) 

r 

° F 


573.0 


572.6 

i 


0.07 

[ - - 


557.4 


557.1 


0.054 


T (MAX/ 

MIN) 
° F 


499 . 8 


i 

500.5 


0.14 


509.2 


510.1 


0.178 


t av ( ss) 

°F 


500 . 4 


501.2 


0.16 

l _ . 


508.2 

i 


509.0 


0.156 


T (MAX/ 

MIN) 
° F 


465 . 3 


464.9 


0.086 


: 486.2 


486.0 


j 0.042 


T s (SS) 
°F 


466.0 


465.5 

3 


0.107 

L ... - \ 


\ 485.2 


484.9 


0.0625 


SETTLING 

TIME 


51.5 


47.4 


8.6 


| 49.5 


52.0 


5 


* 


1.199 


1.2 


0.083 


0.801 


0.8 


0.013 



TABLE (4-1) . COMPARISON OF ANALOG AND DIGITAL SOLUTION 
FOR A STEP LOAD CHANGE OF 20 PERCENT. 
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A i p = +0.4 


A ij; = -0.4 


ANALOG* 


DIGITAL 


% ERROR 


ANALOG* 


DIGITAL 


% ERROR 


K (MAX/MIN) 
$ 


0.117 


0.119 


1.68 


-0.188 


- 0.198 


5.0 


N (MAX/MIN) 


1.40 


1.38 


1.45 


0.586 


0.563 


4.6 


N (SS) 


1.38 


1.37 


0.73 


0.618 


0.595 


4.0 


T (MAX/ 

MIN) 

°F 


580.1 


580.1 


0.0 


548.6 


548.0 


0.123 


T p (SS) 

°F 


580.1 


579.9 


0.034 


549.4 


548.8 


0.123 


T (MAX/ 
AV MIN) 

°F 


495.7 


496.3 


0.12 


515.0 


516.2 


0.232 

i 


t av< ss) 

° F 


497.0 


497.5 


0.10 


512.2 


513.0 


0.156 


T (MAX/ 
b MIN) 

°F 


456.0 


455 . 3 


0.154 


498.0 


495.0 


0.60 


T s (SS) 
°F 


457.0 


456.4 


0.133 


495.2 


495.2 


0.0 


SETTLING 

TIME 


54.2 


50.0 


8.4 


108.0 


98.0 


10.0 


'P 


1.397 


1.4 


0.215 


0.604 


0.6 


0.667 



TABLE (4-2) . COMPARISON OF ANALOG AND DIGITAL SOLUTION FOR 
A STEP LOAD CHANGE OF 40 PERCENT. 
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V. CONTROL 



The purpose of an external control scheme for a nuclear 
power plant must be to set up the desired steady state 
operating conditions, and to restrict transients in vital 
plant variables to be consistant with design limitations. 

The above purposes should be obtained without any significant 
reduction in the inherent stability characteristics of the 
plant without external control. 

Examination of the basic plant transient responses in 
Appendix B indicates that the system under study is very 
stable when subjected to load changes. It is however, 
desirable, as pointed out below, to limit the excursion of 
certain plant variables from their steady state normal value. 
In a plant with a large degree of inherent stability direct 
control of plant variables may be attempted. If a plant has 
poor inherent stability, neutron-level control should also be 
used [ 3] . 

A. STEADY STATE PROGRAMMING 

The pattern that the plant temperatures, pressures, and 
flow rates assume as a function of power-demand is oftened 
referred to as a program [4] . A specific control strategy 
can be devised to obtain the desired program. The two most 
significant programs for discussion are the constant T^ v and 
constant T g programs. 

/ 
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I 



1. Constant T 




Fig. (5-1) . SKETCH OF TEMPERATURE AND PRESSURE 
PROFILES FOR A CONSTANT T^ y PROGRAM. 



Figure (5-1) is a graphical representation of the 
constant T^ v program. At zero power demand, all significant 
plant temperatures (except T ) are at the same value. As the 

r 

power demand increases, the steam temperature and hence 
pressure decrease, while T u and T^ vary in such a manner as 

n L 

to keep T av constant. 




Fig. (5-2) . SKETCH OF TEMPERATURE AND PRESSURE 
PROFILES FOR A CONSTANT T g PROGRAM. 
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Again at no load all pertinent temperatures are equal. 
As the load increases the temperature and the pressures follow 
the pattern as shown in Figure (5-2) . 

3. Comparison of Constant T^ v and Constant T g Programs 

The constant program is the one preferred by the 

reactor. The principle advantage of the constant T^ v program 
for this type reactor is that the primary coolant pressure 
also remains constant. Therefore the need for an external 
pressurizer in the primary circuit is reduced. 

One additional advantage occurs if the plant under 
study had been fuel with an uranium alloy vice an oxide of 
uranium. For a metal alloy fuel the feedback of internal 
reactivity associated with changes in fuel temperature is 
insignificant when compared to the feedback associated with 
changes in average coolant and therefore can be ignored [3]. 
For this case equation (2-9) becomes: 



K 



K R + 



a T,. tT 
c AV 



(5-1) 



Since in stead stake K must equal zero, it is apparent 
from equation (5-1) that there is no need for control rod 
motion to maintain T^ v constant. However the plant under study 
is fueled with UC^ and changes in fuel temperature have a 
strong affect on net reactivity. Examination of equation 
(2-9) now indicates that there is a requirement to change 
control rod position to keep the average temperature constant. 

The advantages of the constant T g Program are apparent, 
from Figure (5-2) , in the design of the secondary loop of the 
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plant. A constant Tg program also provides a constant steam 
pressure which allows for an optimum design of the steam 
plant with conventional controlling devices. 

The two primary disadvantages of the constant T g 
program are first, allowing T^ v to vary over a large range 
(required to keep Tg constant) may cause serious problems in 
the primary coolant pressure system and second, when the plant 
is changing from one steady state power level to another, the 
plant must fight the tendency of the coolant negative temper- 
ature coefficient to hold T, tt constant. 

AV 

With the exception of Nautilus, which uses a constant 
T c program, all Naval plants use a constant T^ v program [16] . 

The N. S. Savannah utilizes a control plan for a constant T... 
program with a plus or minus 3°F dead zone about the normal 
load T av value [17] . 

B. DESIGN OF THE EXTERNAL CONTROL SYSTEM 

Consistant with the advantages of the constant T^y program, 
the simulated power plant will be controlled under a constant 
T^y program. As noted in section (5 A- 3) control rod movement 
is required in order to obtain the desired program. The next 
decision to be made is at what rate may reactivity be externally 
introduced into the reactor. The answer to this and other 
design questions are discussed below. 

1 . Reactivity Rate 

The one basic underlying consideration which dictates 
the design of a control system for a nuclear power plant is 
safety. As shown in Figure (3-3) the power plant operating 
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at normal power can withstand a step reactivity input of 
magnitude 6 k = .469$. This however is not the case for a 
reactor operating at reduced power levels where the negative 
reactivity feedback effect is reduced. Therefore in order to 
provide a maximum safety margin to avoid startup or low power 
accidents input reactivity rates will be restricted to a value 
which can be handled by a reactor operating at lower power. 

The solution to equations (2-7) and (2-8) for ramp input of 
reactivity, K(t), does not give a valid solution to the 
problem because there is no feedback considered. Externally 
introduced reactivity rates are determined largely from 
empirical data on the individual plants involved. A reasonable 
reactivity rate for a plant of this size is lOC/sec (0.1 $/sec) 
[7] . 

It is feasible that this rate of reactivity may be 
controlled by a manual means. However, it was desired that 
the system be fully automated therefore an automatic controller 
was used. 

2 . Automatic Controller 

The automatic controller consists of a hybrid combina- 
tion of an analog control rod positioning system and a digital 
system for comparison of T^ v plus the time derivative of T^ v 
with the reference value of 505°F. It was decided to use the 
+ 3°F dead band error detector of N. S. Savannah as a valid 
method to obtain a discontinuous controller. The principle 
advantage of a discontinuous controller is that it utilizes 
the reactor inherent stability characteristics to reduce wear 
on the control rod mechanisms. The dead zone, although it 
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must be consistent with, the design temperature (i.e., pressure) 
fluctuations in th.e primary loop, discriminates against noise 
on the T av signal. A block, diagram description of the entire 
simulation, plant plus controller, is represented in Figure 
(5-3) . The design of the control rod positioning system and 
a discussion of the hybrid simulation along with transient 
response curves for a step load input of +40 percent of the 
system with external control are shown in Appendix D. 

3 . Evaluation of External Controller 

The external controller was implemented and incorpor- 
ated into the simulation. The system was subjected to step 
load disturbances of plus and minus 20 and 40 percent. The 
improvement in system performance can be seen from tables (5-1) 
and (5-2) . The reduction of the excursions of T^., both 
transient and steady state, and the reduction of settling time 
were considered significant. 

The plant was subjected to step load changes of up to 
80 percent (an unrealistic change) . It was found that the 
system was stable and that the system response improved up 
to step load inputs of 60 percent. Above 60 percent, the 
system response deteriorated and eventually went into a limit 
cycle when the step load disturbance became greater than 75 
percent. The deteriorated performance and eventual instability 
of the system were attributed to the fact that the derivative 
of T^ v became very large and began to work in the opposite 
direction of the T^ v excursions rather than in normal anti- 
cipatory mode of a position plus derivative feedback controller. 
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POWER PLANT 




52 



Fig. (5-3) . BLOCK DIAGRAM OF POWER PLANT WITH EXTERNAL 
CONTROLLER. 







A ip = +0.2 


A ip = -0.2 


EXTERNAL 

CONTROLLER 


NO 

EXTERNAL 

CONTROLLER 


EXTERNAL 

CONTROLLEr 


NO 

EXTERNAL 

CONTROLLER 


* t av 


- 3 ° F 


1 

Ln 

• 

K> 

o 


+ 4 . 2 ° F 


+ 4 . 2 °F 


T av (SS) 


505 .6°F 


500 . 4°F 


505 ' °F 


508. 2°F 


K (MAX/MIN) 


0 .12$ 


0.068 

J 


- 0.16$ 


- 0.081 


PS 

<1 


+ 0.1$ 




o 

• 

o 




SETTLING 

TIME 


36 sec 


51.5 sec 


40.0 sec 


49.5 sec 



TABLE (5-1) . COMPARISON OF SIMULATED POWER PLANT 
WITH AND WITHOUT EXTERNAL CONTROLLER 
FOR A STEP LOAD CHANGE OF 20 PERCENT. 
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A rp = +0.4 


A ip = -0.4 


EXTERNAL 

CONTROLLER 

L. 


NO 

EXTERNAL 

CONTROLLER 


EXTERNAL 

CONTROLLER 


NO 

EXTERNAL 

CONTROLLEF 


AT 

AV 


- 5.1 


- 9 . 3° F 


+ 7 . 2°F 


+10 .0°F 


T av <SS) 


508. 0°F 


497. 0°F 


507. 4°F 


512. 2°F 


K ( MAX/MI N) 


0.17$ 


0.117$ 


- 0.32 


- 0.188$ 


ak r 


+ 0.6$ 




+ 0.6$ 




SETTLING 

TIME 


36 sec 


54.2 sec 


77 sec 


108 sec 



TABLE (5-2) . COMPARISON OF SIMULATED POWER PLANT 
WITH AND WITHOUT EXTERNAL CONTROLLER 
FOR A STEP LOAD CHANGE OF 40 PERCENT. 
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The controller was considered adequate for all 



reasonable power changes . 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



The following conclusions and recommendations are offered 
as a result of the investigation contained in this report: 

A. CONCLUSIONS 

1. A properly designed pressurized water nuclear reactor 
(i.e., one with a negative temperature coefficient of 
reactivity) possesses a great deal of inherent stability when 
the associated power plant is operating in the normal power 
range. It was observed that the negative temperature coeffi- 
cient of reactivity served to stabilize the reactor for dis- 
turbances in both load and input reactivity. 

2. The naturally occuring phenomenon of delayed neutrons 
greatly simplifies the control problems associated with 
nuclear reactors. The effect of this phenomenon is more 
evident at the lower power levels before the temperature 
reactivity feedback effect has become affective. 

3. The control of a stable nuclear power plant operating 
in the power range under a constant T^ v program with constant 
coolant flow rate is a relatively simple task. If sufficient 
stability exists in the reactor then direct control of a plant 
parameter (T^) can be attempted. 

4. If step load changes greater than 60 percent, for which 
the above control system is not suitable, are considered 
feasible a comparator could be installed in the differentiating 
circuit to disregard the derivative if it becomes too large. 
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B. RECOMMENDATIONS 



Further workers in this area should consider the following 
facets of reactor systems. 

1. The above study was concerned with the reactor power 
plant dynamics when excited by both load disturbances and 
changes in reactivity investigated separately. Another area 
of interest is the response of the reactor under the condition 
of load changes and externally induced reactivity changes 
(both aiding and opposing) occuring simultaneously. This 
situation might represent an accidental condition. 

2. The T^ v program used to control the power plant does 

T h + T 

not consider the maximum value to T„. Since T„ tt = x — — 

H AV 2 

the value of T^ might become too high for design considera- 
tions. To institute a control program with a constant T^ v 
and also limiting T^ to some maximum value would make the 
requirement of a variable coolant flow rate necessary. This 
added requirement would necessitate a new mathematical model. 
Under the assumption of a constant flow rate, the time coolant 
spends in the reactor (Tq) is a constant. However, in the 
model for a plant where a maximum value of T^ must be 
considered t» would be a function of T . 

U n 

3. More sophisticated control problems, than reactor 
control at operating conditions, exist in the nuclear power 
plant control field. These are the problems associated with 
the reactor in the startup range and problems involved in the 
mechanical control of the various plant components of the 
secondary loop. The problems in the control of the secondary 
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loop result from the fact that when the reactor is controlled 
under a constant program, the steam temperature and 

pressure vary significantly over the operating range of the 
plant. It is recommended that further investigators consider 
the control problems of the secondary loop. 
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APPENDIX A 



SPECIFIC SYSTEM EQUATIONS 

In the development of the generalized system equation, 
equations (2-28) and (2-34) are expressed in terms of the 
variable T fi . T fi is not a state variable of the system. The 
equation for the primary side of the heat exchanger (2-31) can 
be solved for T £ and this value substiuted into the equations 
(2-28) and (2-34) to place them in a more suitable form. 

This yields : 




The system parameters for the hypothetical plant are stated 
below. It should be emphasized that these parameters do not 
represent any actual real system but that the parameters used 
here are obtained from various sources in the literature [5, 
6, 7, 8, 9] and represent typical values which might be 
associated with a plant the size of the power plant of N. S. 
Savannah . 

1) Fraction of delayed neutrons 8 = 0.0064 

2) Average delay constant of the 6 group delayed 
neutrons X = 0 . 1 
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-3 

3) Mean neutron life time % = 10 sec 

4) Steady state temperature 



A) 


Average fuel temperature 


i-3 

O 

II 


565°F 


B) 


Average coolant temperature 


1 AV C°> = 


505°F 


C) 


Boiler inlet temperature 


t bi (°) = 


520 °F 


D) 


Steam temperature 


T s <°) ‘ 


475°F 


E) 


Cold leg temperature 


n 

o 

ii 


490° F 


Temperature coefficient of reactivity for 


the fuel 


(UO 


2 ) 






a F 


= -1.57 x 10~ 2 $/°F 







6) Temperature coefficient of reactivity for the coolant 
(pressurized water) 

a c = -3.14 x 10~ 2 $/°F 



7) Time constants for the various heat transfer equations 

A) t 1 = 2.0 sec. 

B) t 2 = 0.1 sec. 

C) x ^ = 2.0 sec.* 

D) = 2.0 sec.* 

E) T,- = 5.0 sec.* 

8) In the region of interest the steam pressure-steam 
temperature relation is approximately constant 
2Po 

= 5 psi/°F 



2T, 



9) The values of K R (0) , A^ , K T , were determined 

by setting the equation in which these terms appear to 
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* 






zero (steady state) and solving for the unknown 
variables. The results are: 



A) 


V°> 




24.7275 $ 


B) 


A i 


= 


30.0 


C) 


E 


= 


2.0 


D) 


k t 


= 


1.33 x 10 


E) 


T o 


= 


0.05 sec. 


F) 


T 2 


= 


0.02 sec 



* x 3' t 4' an< ^ t 5 were obtained by engineering approxi- 
mation and are substantially less than those mentioned 
for a stationary power plant [9] . 
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APPENDIX B 



TIME RESPONSE OF POWER PLANT 

As a representative example of system time response 
Figures (B-l) , (B-la) , (B-2) and (B-2a) of this appendix show 

the time response of selected system variables to a step load 
change of +40 percent normal load. Figures (B-l) and (B-la) 
are solutions using the digital simulation. Figures (B-2) 
and (B-2a) are solutions utilizing the analog simulation. 
Pertinent values of the variables are shown adjacent to their 
time response. The system response for other experimental 
load changes are similar in nature to those shown. 
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K (MAX) 



.119 $ 



40 80 120 

NET REACTIVITY VERSUS TIME 



160 secs 



/. 




N (MAX) 
N (SS) 



1.38 

1.37 



40 80 120 

NEUTRON POPULATION VERSUS TIME 



160 secs 



lo 



Ip = 1.0+0 .4y (t-20) 



40 

LOAD CHANGE 



80 



120 



160 secs 



Fig. (B-l) 



TIME RESPONSE OF SIMULATED POWER 
PLANT TO A 40 PERCENT STEP LOAD CHANGE 
(DIGITAL) . 
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T g (MIN) 



455 .3°F 



V7i*V 




0 40 80 120 

STEAM TEMPERATURE VERSUS TIME 



S Wf 



0 40 80 120 

FUEL TEMPERATURE VERSUS TIME 



SOf'F 



160 secs 



T p (MAX) 



T f (SS) 



160 secs 



t av (min) 



t av ( ss) 



0 40 80 120 160 secs 

AVERAGE COOLANT TEMPERATURE VERSUS TIME 



Fig. (B-la) . TIME RESPONSE OF SIMULATED POWER 
TO A 40 PERCENT STEP LOAD CHANGE 



456. 4°F 



580 .1°F 
579. 9°F 



= 496. 3°F 
= 497. 5°F 



PLANT 
(DIGITAL) . 
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K(MAX) = .117 $ 
o 



net reactivity versus time 



N(MAX) =1.40 
1,0 N(SS) = 1.38 





neutron population versus time 
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Fig. (B-2). TIME RESPONSE OF SIMULATED POWER PLANT 

TO A 40 PERCENT STEP LOAD CHANGE (ANALOG) . 
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APPENDIX C 



PREPARATION OF SYSTEM EQUATIONS 
FOR USE ON ANALOG COMPUTER 

Equations (4-2) and (2-36) through (2-42) describe the 
power plant in terms of system variables. These system vari- 
ables have units such as degrees Fahrenheit and dollar units 
of reactivity which must be related to analog computer vari- 
ables which are voltages . In addition to the conversion of 
units, the analog voltages have a maximum value of 100 volts 
so scaling of the variables may also be required. 

To accomplish this conversion a scaling factor coefficient 
x is introduced. The relation between computer variables and 
system variables then becomes : 



N = x n N 
D = x d D 




and so forth. 

The individual x's are equal to the maximum expected value 
of the system variable divided by the maximum desired voltage 
for the computer variables. For example: 

T F = X TF T F 
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100 volts 



assuming T p (MAX) = 600°F and T p (MAX) = 

T 



X TF = *p = 6°F/volt 
F 



A summary of the scaling coefficients is as follows 



x N = . 1 volt 



-1 



X D = . 1 volt 



-1 



x R = .2$/volt 



X KR = * 7 $/ volt 



X = . ] volt 



-1 



X TF X TS X TAV X TBI 6°F/volt 

The system equations are now altered by replacing the system 
variables with their corresponding computer variables and 
conversion factors. The initial condition of the differential 
equations are determined as follows: 

_ T (0) 

V 0) = -i — 

TF 



In a final form the equations now become: 
1. Reactor Kinetics 



N = 



5D 

5-K 



D = 0.1 (N-D) 



K = 3.5 K r - 0.471 T p - 0.942 T^ y 
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Where 



N(0) = 10 volts 
5(0) = 10 volts 
K r ( 0) = 35.36 volts 

2. Heat Transfer from Fuel to Coolant 
T F = - 0.5 (T f - T av ) 5.0 N 

Where T F (0) = 9 4.17 volts 

3 . Heat Balance for Coolant 



'AV 



= -50 



T + T + 
AV F 



40 T, 



Where T Ay (0) = 84.17 



4. Time Delays 



T BI 0,5 ( ~ T BI + 2T AV t c ) 



T c = 0.5 (-T c + 2/3 T g + 1/3 T B][ ) 



Where T BI (0) = 86.67 volts 
T^(0) = 81.67 volts 

5 . Heat Balance for Steam Generator Secondary-Side 
T g = (0.133 T^ - .1324 T^ - 0.0027 T^) 

Where Tg(0) = 79.17 volts 
ip CO ) = 5 volts 
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APPENDIX D 



CONTROLLER SIMULATION 

Analytical design of the external controller could have 
been approached in one of two ways. First a small signal 
linearization of the plant dynamics could have been attempted. 
This method was not valid since two of the plant variables, 

N and D, change in almost direct proportion to load changes. 

For example a 40 percent load change causes a 40. percent 
change in N and D. This change in magnitude of the variables 
was outside the range of small signal application. Second a 
describing function could have been obtained utilizing sinu- 
sodial inputs [18] . The describing function technique is a 
laborious method of obtaining insight into the absolute 
stability limits of the system but gives little insight into 
the controller design problem. In this initial study the 
controller design and evaluation was carried out through 
simulation. If a more detailed study had been undertaken, a 
describing function would have been necessary to verify the 
simulation results. 

The simulation of the external controller was accomplished 
by use of a trial and error technique utilizing the hybrid 
computing feature of the XDS-9300 and Comcor CI-5000 computers. 

The analog voltage representing T^ v was sampled at a 
frequency of 10 samples per second (this sampling frequency is 
well in excess of that required by the sampling theorem) . The 
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values of T were stored in an array. When the above array 
had accumulated 10 members, an approximation to the time 
derivation of T AV (10) - T AV (1) (divided by unit time) . The 
inclusion of the derivative feedback will give the system an 
anticipatory response. After the subtraction the array was 
updated so that a new derivative was available every 0.1 
second. Due to the fact that was changing at a rate of 

less than 0.5 degrees per second at distrubances below 60 
percent, the approximation of the derivative of T ^ was con- 
sidered accurate. The sum of T^ and T^ were subtracted from 
T^y references (505°F) . If the error were more than 3°F an 
error signal of + 1 volt was sent to the analog computer. 



+ V 


K 

in 


+ 0 


WORM GEAR 


± «k r 




S (St +1) 
m 




ASSEMBLY 





Fig. (D-l) . BLOCK DIAGRAM OF CONTROL ROD 
POSITIONING SYSTEM. 



The + 1 volt error signal was placed into a standard 

K 

electrical motor with transfer function — m . t n • As 

S (St + 1) 
n 

indicated in Figure (D-l) in the actual plant this motor would 
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drive at constant speed through a worm gear assembly and 
remove or insert the control rods at a constant speed. The 
gear assembly represents only a linear gain and can be com- 
bined with the motor gain to give a new transfer function 



K. 



lm 



S(ST m + 1) 



The control rod effectiveness (reactivity per 



inch) is a non-linear function of position and is different 

for individual reactors. It is assumed here that the rods 

are operating in a linear mode where (SK^ = .1 $/sec. Under 

the assumption that x = .1 sec. the final transfer function 

m 

becomes g (s + IQ) * Simulation of the drive motor yields the 



transient responses of Figures (D-2) and (D-3) . From the 
simulation of the control rod positioning device the rise 
time is 0.25 seconds and the steady state SK is 10.15C/sec. 

K 

The rise time could be shortened and the steady state K 

K 

remain constant if the forward gain were increased and a 
tachometer feedback loop used. However this would allow the 
possibility of an excessive 6 k r if the tachometer loop were to 
fail. Since the rise time of the uncompensated system was 
quite satisfactory for this application the tachometer feedback 
loop was not employed. 

Transient response of the system with external controller 
is shown in Figure (D-4) and (D-4a) . 
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NET REACTIVITY VERSUS TIME 




N (MAX) 
N (SS) 



1.55 



1.4 




A ip 



+ 0.4 



LOAD CHANGE 

Fig. (D-4 ) . ’ TIME RESPONSE OF POWER PLANT WITH 
EXTERNAL CONTROLLER. 
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T g (MIN) = 
Tg(SS) = 



STEAM TEMPERATURE VERSUS TIME 




T p (MAX) = 
T p (SS ) = 



FUEL TEMPERATURE VERSUS TIME 




t av (min) 

t av ( ss) 



AVERAGE COOLANT TEMPERATURE VERSUS TIME 

Fig. (D-4a) . TIME RESPONSE OF POWER PLANT WITH 
EXTERNAL CONTROLLER. 
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465. 4°F 
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= 508. 0°F 



76 



BIBLIOGRAPHY 



1 . Harrier, J. M. , Nuclear Reactor Control Engineering , 

D. Van Nostrand Company, Inc., 1963. 

2. Glasstone, S. And Edlund, M. C., The Elements of Nuclear 

Reactor Theory , D. Van Nostrand Company, Inc., 1952. 

3. Schultz, M. A., Control of Nuclear Reactors and Power 

Plants , McGraw-Hill Book Company, Inc. , 19 61. 

4. King, C. D. G. , Nuclear Power Systems , The MacMillian 

Company, 1964. 

5. Landis, J. W. , "The Power for the First Nuclear Merchant 

Ship N. S. Savannah," Journal of American Society of 
Naval Engineers , V. 70, p. 629-641, November 1958. 

6. Paulson, A. E. and Smith, W. R. , Observed and Predicted 

Performance of N. S. Savannah Power Core I , Paper 
presented at the American Nuclear Society Topical 
Meeting on Nuclear. Performance of Power-Reactor Cores , 
September 26-27, 1963. 

7. Babcock and Wilcox Company Report 1203 N. L. , Extended 

Zero Power Test- N. S. Savannah Core I , MacKinney, 

A. L. (ed.), January 1961. 

8. The Babcock and Wilcox Company Report 1243, Consolidated 

Nuclear Steam Generator for Merchant Ship Application- 
A Conceptual Design , August 1962. 

9. Duncombe, E. and Rathbone, D. E., "Optimization of the 

Response of a Nuclear Reactor Plant to Changes in 
Demand," IEEE Transactions on Automatic Control, 
p. 277-282, June 1969. 

10. Cameron, W. D. and Taylor, G. R. , "Examples of Analog 

Computer Applications in Nuclear Processes," Simulation , 
p. 23-38, April 1964. 

11. Soodak, H. and Campbell, E. E., Elementary Pile Theory , 

Wiley, 1950. 

12. Doolittle, J. S. and Zerbaw, A. H. , Engineering Thermo- 

dynamics , International Textbook Company, 1964. 

13. Brookhaven National Laboratory Report 815 (T-315) , Non- 

linear Kinetics of a Two Temperature Reactor, by 
L. M. Shotkin , August 1963. 



77 



14. Birch, A., The Kinetics of Pressurized Water and Other 

Reactor Power Plants , Paper presented at the Second 
United Nations International Conference on the Peace- 
ful Uses of Atomic Energy, Geneva, Switzerland, 
September 1-13, 1958. 

15. Weaver, L. E., Reactor Dynamics and Control , p. 50, 

American Elsevier Publishing Company, Inc., 1968. 

16. Sterling, A. G. , "Basic Control Concepts for Pressurized 

Water Reactor," Journal of American Society of Naval 
Engineers, V. 72, p. 693-697, November 1960. 

17. Crouch, H. F. , Nuclear Ship Propulsion , Cornell Maritime 

18. Thaler, G. J. and Pastel, M. P. M., Analysis and Design 

of Nonlinear Feedback Control Systems, p. 139-227, 
McGraw-Hill, 1962. 

19. Coughren, L. D. and Others, "Analog Simulation of Hanford 

Dual-purpose Reactor Plant," Simulation , p. 393-404, 
December 1965. 

20. Harrer, J. M. and Deshong, J. A. Jr., "Discontinuous 

Servo for Control of Power Reactors," Nucleonics , 

V. 12, p. 44-51, January 1954. 

21. Rajagopal, V. and Gallagher, J. M. Jr., "Some Applications 

of Dynamic Measurements in Pressurized Water Reactor 
Nuclear Power Plants," IEEE Transactions on Nuclear 
Science , p. 15-21, April 1967. 

22. "Savannah Control Rod Operation," Nuclear Engineering , 

V. 5, p. 447-454, October 1960. 



78 



INITIAL DISTRIBUTION LIST 



No. of Copies 



1. Defense Documentation Center 2 

Cameron Station 

Alexandria, Virginia 22314 

2. Library, Code 0212 2 

Naval Postgraduate School 

Monterey, California 93940 

3. Professor A. Gerba 4 

Department of Electrical Engineering 

Naval Postgraduate School 
Monterey, California 93940 

4. LT H. Ward Clark, Jr., USN 1 

31 Upper Mountain Avenue 

Montclair, New Jersey, 07042 

5. Dr. G. J. Thaler, Code 52Tr 1 

Department of Electrical Engineering 

Naval Postgraduate School 
Monterey, California 939340 

6. Mr. R. L. Limes, Code 52Ec 1 

Department of Electrical Engineering 

Naval Postgraduate School 
Monterey, California 93940 



79 



Security Classification 



DOCUMENT CONTROL DATA - R & D 

fSrcurify c/e s si tic at ion of title, body of abstract and indexing annotation must be entered when the overall report is classified) 


i OR’C'NA ting activity (Corporate author) 

Naval Postgraduate School 
Monterey, California 93940 


2«. REPORT SECURITY CLASSIFICATION 

Unclassified 


2b. GROUP 


3 REPORT TITLE 

Study of a Marine Nuclear Propulsion Plant Through Computer 
Simulation 


4 DESCRIPTIVE NOTES (Type of report andjnclusive dates) 

Masters' Thesis. June 1971 _ 


5 AUTHORiSi (firsl name, mtddfe initial, la st name) 

Hiram Ward Clark, Jr. 


A REPOR T O A TE 

June 1971 


la, TOTAL NO. OF PAGES lb. NO OF REFS 

81 22 


a a contract or grant no 
b. PROJEC T NO. 

C. 

d. 


9a. ORIGIN A TOR’S REPORT NUMBER!*) 

9 


9b. OTHER REPORT NOSSJ (Any other numbers that may be assigned 
this report) 


10 DISTRIBUTION STATEMENT 

Approved for public release; distribution unlimited. 

2 , 


11 SUPPLEMENT AR Y NOTES 


12 SPONSO R IN G Ml LI T AR Y ACTIVITY 

Naval Postgraduate School 
Monterey, California 93940 



13 abstract 



The dynamics of a proposed marine nuclear propulsion plant, under 
the condition of sudden changes in load, are studied by both an 
analog and digital simulation. Two inherent reactor characteristics, 
internal reactivity feedback and the existence of delayed neutron 
emitters, which effect power plant stability and controllability are 
considered in detail. 

Comparison of the simulation results with practical data indicates 
that the simulation represented the dynamics of a marine nuclear power 
plant with a sufficient degree of accuracy to be extremely useful in 
the study of marine nuclear power plants. 

Evaluation of the simulation results indicate that the proposed 
power plant is inherently stable when operating in the normal region. 
The desirability of control schemes are discussed and a control scheme 
utilizing a constant average coolant temperature program is imple- 
mented. This external control scheme significantly improved the 
response of the system. 



DDtr.,1473 IP#GE " 



S/N 0101 -807-681 1 



80 



Security Classification 



A-31408 




Security Classification 



key words 



Marine Nuclear Power Plant 
Computer Simulation 



ROLE W T 



, F °o R v M J473 (back, 

101 * 607-6921 



81 



Security Classification 



A - 3 | 409 



18 APR 73 



I 




Thes i s 

C4818 

c.l 






^D'ltlDERY 
2 0 7 0 1 



1 

Thesis 132941 

C4818 Clark 

c.l Study of a marine 

nuclear propulsion 
plant through com- 
puter simulation. 

* D'l HDERY 

19 APR 73 2 0 7 0 1 



132941 

Clark 

Study of a mari ne 
nuclear propulsion 
plant through com- 
puter simulation. 



thesC48l8 

Study of a marine nuclear propulsion pla 




I II I II II I ll i ii i 

3 2768 001 02602 4 

DUDLEY KNOX LIBRARY 













